%--------------------------------------------------------------------------
% Script to overlay percentiles
%--------------------------------------------------------------------------

clear; 
clc;
close all

set(0,'defaultTextInterpreter','latex');

% set specs 
xmin = 0;
xmax = 2.6;
xn = 100;
xgrid = linspace(xmin, xmax, xn);

nfVARSpec = '1';
nKSpec    = 'K8_';
sName = ['fVAR', nfVARSpec];

% load sample percentiles (from the data)
dataDir = [pwd, '/', 'Data/Para1/' ];
sample_percs = csvread( [dataDir, 'percentiles_data.csv'], 1, 1); % sample percentiles

% load estimated percentiles 
estDir = [pwd, '/', 'Results/Para1/' sName, '/'];
estimated_percs = csvread( [estDir, nKSpec, sName, '_PredPctl_MLE.csv'], 1, 0); % estimated percentiles

%--------------------------------------------------------------------------
% Percentiles Comparison 
%--------------------------------------------------------------------------
T = length(sample_percs(:,1));

figsaveDir = [pwd, '/', 'Figures' ,'/', 'Para1', sName,'/'];
[~, ~, ~] = mkdir(figsaveDir);

%start_period = 1989.25;
%period = linspace(start_period,start_period+0.25*(T-1), T);
period = 1:T;

blue_colors = [0.4 0.4 1; 0 0 0.9; 0 0 0.8; 0 0 0.9; 0.4 0.4 1];
red_colors  = [1 0 0; 0.9 0 0; 0.8 0 0; 0.9 0 0; 1 0 0];

figure(1);clf;
%set(figure(1),'PaperType','usletter','PaperOrientation','Landscape','PaperSize',[11*2 8.5],'PaperPosition',[0.1 0.1 11*2 8.5]);
set(figure(1),'PaperType','usletter','PaperOrientation','Landscape','PaperPosition',[0.1 0.1 11 8.5]);
for ii=1:5
    plot(period, sample_percs(:,ii), 'Color', blue_colors(ii,:), 'LineStyle', '-', 'LineWidth',4)
    hold on
    plot(period, estimated_percs(:,ii), 'Color', red_colors(ii,:), 'LineStyle', '-.', 'LineWidth',4)
    hold on
end
set(gca,'FontSize',30)
xlim([period(1) period(end)])
%ylim([0 3])
%axis([1 T 0 3])
sNameFile = [nKSpec, 'Pctl_MLE_All.pdf'];    
saveas(figure(1), [figsaveDir sNameFile] );
%close all


